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Quantum critical systems derive their finite temperature properties from the 
infiuence of a zero temperature quantum phase transition.^ The paradigm is es¬ 
sential for understanding unconventional high-Tc supercondnctors and the non- 
Fermi liquid properties of heavy fermion compounds. However, the microscopic 
origins of quantum phase transitions in complex materials are often debated. 
Here we demonstrate experimentally, with support from numerical renormaliza¬ 
tion group calculations, a universal crossover from quantum critical non-Fermi 
liquid behavior to distinct Fermi liquid ground states in a highly controllable 
quantum dot device. Our device realizes the non-Fermi liquid two-channel Kondo 
statebased on a spin-1/2 impnrity exchange-coupled equally to two indepen¬ 
dent electronic reservoirs.® Arbitrarily small detuning of the exchange couplings 
resnlts in conventional screening of the spin by the more strongly coupled chan¬ 
nel for energies below a Fermi liquid scale T*. We extract a quadratic dependence 
of T* on gate voltage close to criticality and validate an asymptotically exact de¬ 
scription of the universal crossover between strongly correlated non-Fermi liquid 
and Fermi liquid statesJ30 

A conventional second-order quantum phase transition (QPT) features quantum mechan¬ 
ical fluctuations of a classical order parameter. Some second-order QPTs in heavy fermion 
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materials, notably CeCug-ajAua, and YbRli 2 Si 2 , defy easy description in this scheme, and 
their quantum critical behavior instead appears to be related to the breakdown of Kondo 
screening.® Distinctive non-Fermi liquid behaviors appear above a so-called Fermi liquid (FL) 
scale that vanishes at the quantum critical point (QCP); away from the QCP, a crossover 
from non-FL to FL behavior is observed at low energies. A diverging effective mass m* at 
the QCP, seen in both materials, signihes the absence of quasiparticles at the Fermi surface.® 

In many heavy fermion materials and in high-Tc superconductors, the relevant degrees 
of freedom and the effective Hamiltonian can be controversial. We aim to understand quan¬ 
titatively a second-order QPT outside the usual order-parameter-fluctuation description. 
Quantum dots provide an experimental framework for realizing known quantum impurity 
Hamiltonians that can feature tunable second-order QPTsElE^ However, QCPs are chal¬ 
lenging to reach even in engineered systems, since perturbations that steer away from quan¬ 
tum criticality may be inherently uncontrolled, as in two-impurity Kondo experiments to 

date.™ 

At the QCP of a two-channel Kondo (2CK) system, a single overscreened spin yields a 
non-FL state with no quasiparticles (i.e. only collective excitations) at the Fermi surface. 
An order parameter is typically not invoked; rather, the critical behavior is owing to the 
single spin. A FL scale T* results from several relevant perturbations: Zeeman splitting, 
difference in exchange couplings, and charge transfer between the two channels. Requiring 
that all these perturbations be small would seem to diminish prospects for observing the 
QCP in bulk systems. Nonetheless, two-channel Kondo physics has been invoked to explain 
experiments on heavy fermion material^^^HIS] two-level tunneling centers.^^^®^ A 2CK 
state has been predictecP and observecP in a quantum dot tunnel-coupled to a “metallic 
grain,” an electron reservoir big enough to have a small level spacing A < kT but small 
enough to retain a charging energy Eq ^ kT, at temperatures of interest. The metallic 
grain provides an independent screening channel, as the grain’s charging energy strongly 
suppresses inter-channel charge transfer. Non-FL behavior was observed, as were the FL 
single-channel Kondo states far from the QCP, but the crossover to those FL states was not 
explored. The universal crossover functions were however calculated by NRC.^^ Recently, a 
description of the crossover has been found using Abelian bosonization and conformal held 
theoretical methods, yielding asymptotically exact predictions for conductance in the regime 
where V,T,T* <^Tk^ 

In this work, we show how hue control over the 2CK state in a mesoscopic device allows 
direct comparison to exact results in the crossover regime, yielding T* as a function of gate 
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detuning away from the QCP. The device (Fig. [^) is fabricated by lithographically pat¬ 
terning gate electrodes on a GaAs/Alo.sGao.rAs heterostructure hosting a two-dimensional 
electron gas. The device is abstracted in Fig. [^. Despite the number of gates, the device is 
conceptually simple (Fig. [^): a metallic grain (red) and two leads (blue) are each tunnel- 
coupled to a quantum dot (green) at rates Fc and F, respectively. The charging energy is U 
{Ec) for the dot (grain) (full Hamiltonian in supp. info). In this experiment, two-terminal 
conductance G = dl/dVsd is measured between the pair of leads (Methods sec. [^. We use 
Fbwt (Flp) to tune the grain level 0 (dot level e). 



Figure 1: Device and model, (a) SEM micrograph of a nominally identical device. The five 
brighter features seen coming in from the left are metal bridges suspended above the sample surface, 
(b) Schematic of the device with labeled gate electrodes. Gates BWT, BP, and BWB define the 
grain (red) along with LBT and LBB; the last two also control the dot-grain coupling. Gates LWT, 
LP, and LWB define the dot (green), along with LBT and LBB. Gates BR are used to isolate 
the dot measurement circuit. Other gates are held at a fixed voltage throughout the experiment. 
Gonductance is measured between source and drain leads (blue). The four gray stars indicate 
additional ohmic contacts which are floated during measurement, (c) Model of the system used for 
the NRG calculations. L^ is the dot-grain coupling, L the total dot-lead coupling (sum of couplings 
to source and drain leads). The source and drain leads together act as one channel in the spin 
2CK regime, and the Goulomb-blockaded grain acts as an independent channel. Full Hamiltonian 
in supp. info. 

We first identify the set of QGPs in the {—e/U, —(j)/Ec) plane for hxed F,Pg. For our 
model Hamiltonian, quantum critical “2GK lines” periodic in the grain charge are expected 
instead of isolated QGPsBUHHH! Figure shows the 2GK lines overlaid on numerical renor¬ 
malization group (NRG) calculations of G{—e/U, —cp/Ec) using realistic device parameters. 
We focus on the spin 2GK regime, though charge fluctuations may be important elsewhere.^ 
To directly compare to the experimentally measured conductance data of Fig. |^, Fig. 
adjusts the NRG calculations of Fig. ^ to account for the cross-capacitance between I4 p 
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Figure 2: Quantum phase transitions, a) NRG calculations of G'(14d = 0) for symmetric 
source-drain coupling (T = 20 mK). Parameters: U = 2 meV, P = 0.123 meV, Pc = 0.106 meV, 
Ec = 0.15 meV, bandwidth D = 1 meV. 2CK lines are determined by analysis of the finite 
size spectrum, b) The calculations of a) plotted with an e-dependent shift in cj) and rescaled by 
a constant factor for comparison with c), to account for unequal source-drain couplings. White 
lines indicate 2CK lines, c) Experimentally measured G'(14d = 0). Gates Fbwt and Flp act 
approximately like —$ and —e. The dashed line indicates the cut direction of Fig. [^. d) NRG 
calculations of the equilibrium spectral functions A{u, T) for e, <1> as marked in b). The black trace 
is the spectral function ^ 2 Ck(W)F, hp) from GET {5p = —0.0297r, Tk = 19 ^eV). e) Measured 
G{Vsd,T) for Flp, Fbwt as marked in c). The black trace is Y 2 ck{^/T, 5p)/, rescaled based 
on an estimate of source-drain coupling asymmetry {5p = —O.OIGtt, Tk = 50 //eV). The range in 
(eVsd/kT) decreases as temperature increases because we measure a fixed range in 14 d- 
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and the grain. 

To identify transport signatnres of qnantnm criticality along the 2CK line, we look for 
the characteristic sqnare-root scaling of G{Vsd,T) derived from the CFT of Affleck and 
Lndwig.l^ The CFT yields temperatnre-dependent spectral fnnctions A 2 ck{^-,T, 6 p), where 
6 p is a phase shift from potential scattering. These are closely related to G{Vsd,T) for 
u —)■ —eVsd (Methods sec. [^. A scaling collapse of G{Vsd,T) is expected: 


G{0,T) - GiV,d:T) 


oc —F2CK(-e'Fsd/fcT,hp) 


( 1 ) 


VT 

where Tk is a scale below which the 2CK physics is observed and T 2 ck(— el4d/fcT, 5p) a 
nniversal fnnction closely related to A 2 ck(i^, T, 5p) (Methods sec. |^. 

Fignrej^ shows spectral fnnctions A(a;, T) calcnlated by NRG. Importantly, the spectral 
fnnctions collapse onto A 2 Ck(i^, T,5p), with the horizontal axis scaled to emphasize the 
behavior for large uj/T. Measnred G{Vsd,T) on or very near the 2CK line (Fig. [^) collapse 
similarly, except for the 20 and 40 mK traces at positive Kd- This deviation could result 
from a small T* < Tg, the base electron temperature. Data taken at more negative Vlp show 


very clear 2CK scaling (supplemental info) but are less suitable for analyzing the crossover. 
Experimental Tk ~ 50 peV should only be trusted up to factors of order unity: in Eq. ^Tk 
enters only as a scale factor, and other scale factors like source-drain coupling asymmetry 
must be estimated. 

Having identified the 2CK lines in Fig. we consider how to perturb the quantum 
critical state. In the 2CK model, a single FL scale T* suffices to describe any combination 
of symmetry-breaking perturbations.^ The limit uj,T,T* Tk permits an exact expression 
for the scattering T-matrix in the low-temperature 2CK crossover, found by Sela, Mitchell, 
and Fritzl30 In our experimental conhguration the T-matrix is diagonal: 


27liuTaa,aa(,UJ, T, 6p) = 1 - 6^*^^(2) 

with the universal complex-valued function ^ encoding the crossover physics. These 

diagonal elements relate to A{u,T) and thus to experimental G = dl/dVgd for highly asym¬ 
metric source-drain coupling, u is the bare density of states per spin in the leads, a is the 
spin index, and a = 1 (-1) labels electrons in the leads (grain). The F-matrix gives a (spin 
and channel dependent) scattering phase shift that is a function of the relative strengths of 
any perturbations present. Negligible charge transfer between channels and zero magnetic 
field yields Saa,aa = with -|-(—) indicating the dot is more strongly exchange-coupled 
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to the grain (leads). The factor accounts for additional spin-independent phase shifts 
from potential scattering. We £x Saa,aa = a and let 6 p jump by 7r/2 to account for sign 
changes. 

To observe the FL crossover experimentally, we £x Wp = —236 mV (dashed line in 
Fig. [^) and detune the exchange couplings using Vbwt- Moving slightly away from the QCP 
so that T* ~ Te, we still measure a VT scaling collapse for T > 50 mK (Fig. [^). These 
high T data are £t nicely using the Affleck-Ludwig CFT result with small Sp (black line). 
The clear scaling behavior at high-T can only be observed for VpwT in a small neighborhood 
around the QCP. Below 50 mK, prominent deviations from 2CK scaling develop, which we 
attribute to a crossover into a FL state where the grain screens the dot spin. Near zero bias 
these low-T traces are £t by the crossover theory with similar, small 5p (Fig. [^). We stress 
this is a non-trivial regime since T* ~ Te; asymptotics of the FL hxed point are insufficient to 
describe the observed behavior. For larger |eVd/fcT|^/^, the |eVd|^'^^ dependence of G(Vd )5 
appearing linear on these axes, heralds a return to 2CK behavior. 

Generically, T* should depend quadratically on the strength of symmetry-breaking per¬ 
turbations near the QCP0l^ (Fig. [^). Measured G(Vd,f^wT) reveals periodic zero bias dips 
that transition sharply to zero bias peaks as Vbwt is increased (Fig. |^, top). The zero bias 
dip (peak) corresponds to a T = 0 ground state where the grain (lead) screens the dot spin; 
these are separated by a QCP. In Fig. (middle), T* depends quadratically on Vbwt away 
from the QCP, although the curvature differs between the two sides of the QCP, which have 
different ground states. This quadratic behavior holds over a larger range of Vbwt than 
we might have expected considering that generically the exchange couplings do not depend 
linearly on gate voltage. The phase shift ~ 0 on one side of the QCP, and appears to 
approach 7r/2 on the other (Fig. [^, bottom). Between QCPs, 5p varies smoothly. T* and 5p 
are not plotted directly to the right of each QCP, reflecting the ambiguity of htting a small 
crossover peak on top of the 2CK peak. Both T* and 5p are insensitive to small changes in 
the range of I4d used for htting (supplemental info). 

Many features of these observations are corroborated by htting the crossover theory to 
spectral functions from NRG, which yield conductance via Eq. I (Methods). The NRG 
conductance (Fig. [^, top) shows zero bias dips transitioning into peaks, as well as the shift 
of the peak toward positive —w, as in transport spectroscopy (Fig. [^). The 0-dependence of 
T* (Fig. middle) shows asymmetric parabolas like in the measurements. The extracted 
6 p (Fig. 1^, bottom) reproduces the rapid ti/2 phase shift across a QGP, with an otherwise 
smooth 0-dependence. The 7r/2 shift rehects a sign hip in Saa,aa between distinct FL ground 
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Figure 3: Crossover from quantum criticality to a Fermi liquid. Vlp = —236.0 mV 
for experimental data, a) Measured G{Vsd,T). At Vbwt = —376.4 mV, a thermally broadened 
spectral function from the 2CK CFT (5p = —0.0227r, solid black line) describes the high-T data, 
b) Same data in a). G(14d) at low energies is fit to thermally broadened spectral functions from 
the crossover theory (top: 20 mK, bottom: 40 mK; 5p = —0.0457r, T* = 0.5 /reV). Fitting details 
in Methods, c) Quantum criticality occurs for energies above the Fermi liquid scale T* (gray 
paraboloid), which should depend quadratically on the coupling asymmetry Ji — J 2 between the 
two channels as well as on the Zeeman splitting Ez- We vary T* by tuning Ji — J 2 (cut along red 
parabola), d) Extraction of T* and 5p from measurements. The triangle denotes l^wT for a) and 
b). Top: G(Vsd)F = 20 mK). Middle: T* from crossover theory fits to experimental G{Vsd-,T). 
Red traces are parabolas with T* = 0 at the QCP and unequal scale factors on either side of the 
QCP. The largest T* values may not be much less than Tk, so the crossover theory is not strictly 
valid for all Vbwt- Labels indicate approximate QCP locations. Bottom: 6p from the crossover 
theory fits. Error bars reflect 1 s.d. confidence intervals from the fits, e) Extraction of T* and 5p 
from NRG calculations. Parameters as in Eig. Top: G{—uj), rescaled to match maximum G of 
d). Middle: T*. Bottom: 6p. 
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states, where either the grain or leads screen the dot spin.l^ A perfect correspondence 
between experiment and NRG shonld not be expected, since only U and Ec may be extracted 
directly from measurements. Yet both experiment and NRG are well described by the 
crossover theory, and key experimental features are reproduced in the NRG calculations. 

The experimental and numerical corroboration of analytical results in the vicinity of a 
QGP is a milestone in our understanding of correlated electron systems, with implications for 
high-Tc superconductivity and heavy fermions. An essentially identical universal crossover 
(same Q ^), but different symmetry-breaking perturbations) is expected for the two- 
impurity Kondo model. Future works will address the full phase diagram of the device, 
which may host charge 2GKEI1E11E3 and SU(4) Kondo regimesHlEH Our device geometry 
could enable Aharonov-Bohm interference measurements to probe phase coherence of low- 
lying excitations in the non-FL 2GK statep^l^ giving new insight into the nature of a local 
non-FL. 
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Methods 


1 Measurements 

The measurements are performed in the mixing chamber of a wet dilution refrigerator (Ox¬ 
ford Kelvinox TLM) with a base electron temperature Te = 20 mK, verihed by Coulomb 
blockade thermometry. The device was cooled down with -1-300 mV bias on all gates to 
enhance charge stability by reducing the range of voltage needed for operation. For all mea¬ 
surements we use an SR830 lock-in amplifier with 1 /iV excitation at 33 Hz and a custom 
10® V/A gain current preamplifier (design by H.K. Choi and Y. Chung, see related publi¬ 
cation®^. A custom voltage source with six 20-bit channels and eight 16-bit channels was 
used to control the gate and source-drain bias voltages (design by J. MacArthur, assembled 
and calibrated by AJK). 

The biased source lead in any source-drain bias spectroscopy was determined to be weakly 
coupled to the dot: At zero bias, we pinch off the source lead’s coupling to the dot T^ (e.g. 
using Vlwt) and observe a decrease in the overall conductance scale, without appreciable 
changes in the conductance features after accounting for capacitive shifts from gating. This 
implies that the unbiased drain lead’s coupling to the dot T^ largely determines the total 
dot-lead coupling rate, since T = Ts-t-Td ~ T^, i.e. the dot was nearly in equilibrium. 
In comparing Fig. 2a and Fig. 2c, the ratio of maximum conductances is 0.464. If these 
numerical and experimental data are assumed to be directly comparable, then the asymmetry 
prefactor 4FsFrf/(Fs 4- Tdf = 0.464, yielding F^/F^ = 15%. 

It is well known that applying source-drain bias will cause unintentional gating as a 
secondary effect. This would be deleterious to observing quantum critical behavior, which 
depends sensitively on the dot and grain levels. We compensate for shifts in the grain level by 
compensating changes in 14^ with changes in Vbwt- This compensation can be determined 
easily in the regime e/U > 0 or e/U < —1. We expect the grain level to be much more 
sensitive than the dot level for the same change in energy since Ec <C U. 

2 Fitting range 

When fitting the crossover theory to experimental data, we fit G{Vsd,T) only in a small 
window of I4d of ±6 fiV around zero, regardless of temperature. A priori, T* is unknown 
and it only makes sense to £t I4d < a few T*. Additionally, thermal broadening of high 
energy features can in principle spoil the scaling of the low energy features, even for oth¬ 
erwise sensible ranges of I4d- At minimum the 20 and 40 mK traces are used for fitting. 


11 


but sometimes also the 52 mK and possibly the 70 mK traces, provided T <T* (the htting 
process is somewhat iterative in this respect). Once the temperatures to be used in htting 
are decided for a given value of f^wT, the htting considers data from all of those tempera¬ 
tures simultaneously. Fitting the crossover theory to NRG calculations is done analogously 
(window of u of ±6 /reV about zero). 


3 Relationship of G(V"sd, T) to spectral functions 

The diherential conductance G = dl/dVsd measured from source to drain lead through the 
small dot is a function of source-drain bias and temperature and can be compared directly 
to NRG calculations in case of a strongly asymmetrical source-drain coupling. In the case 
of weak coupling to the biased source electrode (F^ -C F^^), the diherential conductance can 
be related to the equilibrium spectral function as 


G(Kd,T) 


2e^ 4F,F^ 

h (F, + F,)2 




(-eKd),T) 

dco 


A{oj,T). 


( 3 ) 


The asymmetry prefactor is a function of the source and drain couplings, F^ and F^;, and is 
assumed to be much smaller than one. Either lead may assume the role of source or drain. 
The derivative of the Fermi-Dirac distribution f{u,T) is convolved with a spectral function 
Aiu, T) from the 2GK or crossover descriptions. The spectral function can be related to the 
T-matrix: 


A{U, T) = -TIU ^ Im [Taa,aa (w, T)] 

(T 



( 4 ) 


where u is the bare density of states in the leads, cr is a spin index, and a is a channel index 
(we £x a = —1 for the source and drain leads). The T-matrix represents the scattering 
between different states induced by the interaction part of the Hamiltonian and can be 
computed numerically exactly by NRG. It is related to the quasiparticle self-energy.®l 


4 Fitting expressions for 2CK 

In equilibrium, the conduction electrons’ scattering T-matrix is proportional to the self¬ 
energy. In case of the quantum dot system considered here, the latter quantity translates to 
the Green’s function of the d-level of the small dot. This allows us to use the exact F-matrix 
at the 2GK hxed poinlP^ and express the equilibrium spectral function of the small dot in 
the limit T* -C cu, T as 
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A{u,T) A2 ck{!^^T,5p) = Im i ( 1 — 


(5) 



- uf/^2Fi{?>/2, 3/2; 1, u) - -u-F‘^{l - 

TT 


where 2 -F"i(a, b; c, z) is the Gauss hypergeometric function, (3 is inverse temperature, and hp is 
the scattering phase shift. We hx the dimensionless parameter A = —0.09 so that the spectral 
function drops to half of its a; = 0 value at a; = Tk in the limit T — 012211221 Equation (|^ 
immediately implies that {A2 ck{0 -,T,6 p) — A 2 ck{oo., T, 5p))\JTk/T is a universal function 
of cn/T, which when convolved with a Fermi function gives the function Y 2 CK{—&Vsd/kT, 6p) 
of equation ([^. We stress that this u/T scaling is a special property of the 2CK hxed point. 
When htting the experimental data, we shall assume an asymmetrical coupling to the leads 
(see Methods sec. [^. 

5 Fitting expressions for crossover 

At frequencies and temperatures far below the two-channel Kondo temperature Tk, we can 
use the crossover form of the T-matrix derived in RefsHHH to express the d-level’s equilibrium 
spectral function. Here we obtain the following expression: 

A(a;, T) ^ App(n;, T, 6 p) = Im ^ (l - (q, t) ) , (6) 

where 5p ~ 0 ((5p ~ 7r/2) in case the dot is coupled more strongly to the grain (leads), and 



Gi^,T) = 




2 ^ 2nT I 


tanh r 1 -I_ 


2T 


■X 


27rf J 


^oo ^ixiolirT 

dx -^—Re 

smha; 


,11 1 


1 — coth X 


(7) 


is a universal function of rescaled energy 0 = oo/T* and temperature T = T/T*. For 
equation ([^ only, F is the gamma function, not a tunnel rate. Again, when fitting to 
experimental data, the spectral function must be thermally broadened (see Methods sec. [^. 
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Figure SI: SEM micrograph of a nominally identical device (5 kV acceleration voltage). The 
device is tilted 40° with respect to normal incidence. 

SI Device 

The 2DEG is 50 nm deep and has an electron density n = 3.3 x 10^^ cm“^ and mobility 
p = 1.2 X 10® cm^/Vs. 

Figure la shows a top-down SEM micrograph of the device. This view is appropriate for 
labeling the gate electrodes and explaining the function of each gate, but the air bridges are 
hard to see. In Fig. SI we show a view of the device at a 40° tilt with respect to normal 
incidence. The hve air bridges clearly rise above the gate electrodes underneath. The device 
is rotated approximately 180° with respect to the orientation of Fig. la. 

As initially fabricated, the bridges did not make good electrical contact to the gates. 
This problem was remedied with an in-situ platinum deposition procedure to be described 
in a forthcoming publication.^ 
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S2 Hamiltonian 


In our numerical calculations the dot-grain system is modeled by the following Hamiltonian: 


H, 


device 


= H, 


dot 


H, 


grain 


H, 


leads 




tunneling? 


( 1 ) 


where 

Hdot = X] ^ did<T + UfUfhi, (2) 

cr 

describes the dot, with e the on-site energy, d\ the creation operator of an electron with spin 
a and ha = S^da representing the occupation number operator. U is the correlation energy 
between the two electrons residing in the dot. The grain is described by 

digrain ^ ^ T (^5 -^o) T 4^i.d^g -^o) • (3) 

p,a 

The creation operator of a spin-a electron with momentum p and energy Sp in the grain is 
denoted by aj,^, Ec is the charging energy of the grain, while —cf) plays the role of a gate 
voltage, hg is the electron number operator of the grain, hg = Yhp a • Opo-Opo-and Nq denotes 
the number of excess electrons in the electrically neutral grain (0 = 0). The non-interacting 
quasiparticles in the leads are described by: 

hfleads akCa.ka^akcT, (4) 

a={U,L} k,a 


where is the creation operator of a spin-a electron with momentum k and energy Sak in 
the upper [a = U) or lower [a = L) lead. The last term in (1) is the tunneling Hamiltonian, 
which is given by 

hftunneling tQ ^ ^ -|- d^^Clpaj ^ ^ ^ ^ (^akadcr T d^^Cako^ ■ (5) 

P,o- a={U,L} k,a 

The tunnel matrix elements between the leads (grain) and the dot are denoted by ta (tc) 
and are assumed to be independent of momentum. The strengths of the couplings are given 
by Tg = TTz/cItcP and To = nualtal"^, respectively, where VaiT^c) is the density of states 
for lead a (grain). In the numerical renormalization group (NRG) calculations the energy 
spectrum of the grain is assumed to be continuous and the densities of states for leads and 
grain are taken to be constant and equal: I'a = = ^ = 1/(2T*), with D = 1 being the 
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band halfwidth used as the energy unit in NRG calculations. 

In Hamiltonian (1) we have neglected the dot-grain capacitive coupling, which can give 
rise to a term of the form, Udgifig — NQ)hd, where fid = + fii- An estimate for Udg 

is extracted experimentally in Section S4, and is believed to play no role for the present 
analysis. 

S3 Summary of the NRG calculations 

S3.1 NRG calculations 

To solve the Hamiltonian (1) we use the numerical renormalization group method.^’^ First, 
we introduce the collective charge operators (bosonic operators) for the grain, 

OO OO 

N= m|m)(m| and |m±l)(m|. (6) 

m=—OO m=—oo 

Strictly speaking, the identity, N = hg, must be fulhlled, but within the NRG approach 
this constraint can be relaxed by treating N as an independent quantity. This is possible 
as the spectral properties of the system are not sensitive to the exact number of conduction 
electrons present in the grain in the limit of inhnitely small level spacing. To extract the 
hnite size spectrum and determine the location of the two-channel Hondo (2GK) lines, how¬ 
ever, a projection to the physical subspace was necessary. In our calculations we took into 
account seven charges in the grain. Using the above charge operators, the grain part of the 
Hamiltonian (3) can be rewritten as 

^^grain = “ ^o)- (7) 

p^cr 

The operators capture the charging transitions of the grain and enter explicitly in the 
tunneling Hamiltonian (5), which now reads 

FTtunneling = • ( 8 ) 

P,(T h (7 

The hrst term in (8) describes the dot-grain tunneling, while the second term accounts for the 
dot-leads coupling. This second term is obtained by performing an orthogonal transforma¬ 
tion® from the two-lead basis to an effective single lead with resultant coupling T = T^, -t-Tf/. 
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The resulting Hamiltonian consists then of two conduction bands coupled to a complex 
impurity composed of the grain [N) and dot degrees of freedom. 

The core of the NRG procedure is the logarithmic discretization of the conduction band 
with discretization parameter A and mapping of the conduction band onto a semi-inhnite 
chain with exponentially decreasing hoppings. The Hamiltonian can then be diagonalized 
in an iterative fashion. In our calculations we used discretization parameter A = 2 and 
kept 4000 states at each iteration. We also exploited the SU{2) symmetry of the total 
spin and two U{1) symmetries for Ni = hd + hcb + hg and N 2 = hg — N, where ficb is 
the electron number operator in the hrst conduction channel (leads coupled to the dot) 
and hg is the electron number operator in the second channel (the grain). We performed the 
full density-matrix numerical renormalization group calculations (fDM-NRG),^’^^® employing 
the Budapest Flexible DM-NRG code,^ to compute the normalized dimensionless spectral 
function, A(a;,T) = ti{Tl + rf/)Ad(ti;, T), where Ad{u,T) is the spectral function for the dj. 
operators that describe the dot level. The linear conductance through the small dot can be 
then determined with the equation 


G 


2e^ ATuTl f, 
h (Tu + TLfJ 




A{u,T), 


(9) 


where f{u,T) is the Fermi-Dirac distribution function. 


S3.2 Shifting of NRG calculations in Fig. 2b 

In Fig. 2b we incorporate a linear e-dependent shift into 0 to obtain agreement between 
NRG calculations and experiment (Fig. 2c). The agreement is obtained by hrst rescaling 
the NRG calculations so that the maximum value of G is the same. The global scaling takes 
into account the source-drain coupling asymmetry, as explained in Methods. Then, the sharp 
features in the cut taken at Vlp = -260 mV are compared with NRG calculations to establish 
Vlp = —260 mV ~ —e/U = 0.55. Finally, another cut for hxed Wp is taken to establish a 
linear relationship between Vlp and —e/U. The two points give a linear relationship between 
Vlp and —e/U. One global offset in —(f)/Ec then suffices to give good agreement everywhere. 
Using this method we hnd -^/Eq = —0/iVc —3-1 (—e/U —1.5), where -^/Eq is the vertical 
axis of Fig. 2b. Physically, the linear dependence of <h on —e/U can be understood as a 
consequence of the indirect capacitive coupling between Vlp and the grain. 
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Figure S2: Measurement of U. Within each Coulomb diamond we label the number of 
electrons on the dot as determined by charge sensing techniques. The intersections of the 
lines indicate t/ ~ 2.9 meV. 

S4 Extracting device parameters 

From measurements in the Coulomb blockade regime we are able to determine the dot 
charging energy U and the grain charging energy Ec- In determining Eq we justify treatment 
of the grain as a continuum. We also determine bounds on any dot-grain charging energy 
Udg neglected in the model. 

The dot charging energy U = 2.9 meV is determined from source-drain bias spectroscopy 
of the dot (Fig. S2). In previous cooldowns U has varied between 1 and 3 meV, perhaps owing 
to how U depends sensitively on the number of electrons in the few electron regime. We use 
U = 2 meV as the model parameter in NRG calculations, and note that the calculations 
should be relatively insensitive when U > D = 1 meV, the electronic half bandwidth used in 
calculations. This value of D corresponds roughly to the internal level spacing on the small 
dot, providing a high energy cut-off. 

The grain charging energy Ec = jC = 160 /reV is measured by source-drain bias 
spectroscopy of the grain (Fig. S3). We compare this measurement to geometric estimates. 
A common rule of thumb is that upon gate depletion, the extent of the depletion region 
extends as far from the gate laterally as the 2DEG is deep. This means that the area of the 
grain should be the area outlined by the gates, less some area in a ~ 50 nm dead region 
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-197 VBWT(mV) -195 


Figure S3: Measurement of Ec- a) Measurement scheme. G(i4d) is measured using the 
grain’s own pair of measurement leads, which are isolated from the measurement leads of 
the dot by depleting gate BR. Gate BL is depleted to avoid shorting conductance through 
the channel just left of the grain, b) G(14d) through the grain in the Coulomb blockade 
regime (T = 20 inK). The intersections of the lines indicate Eq ~ 160 /reV. 
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Figure S4: Bounding Udg from measurements in the Coulomb blockade regime. (a) 
C(14d, Rlp) through the dot in the Coulomb blockade regime, with both the dot and 
grain formed. Here Vbwt is such that the grain is Coulomb blockaded. From the slopes 
of the dashed lines overlaid on the peaks in the data, we determine lever arms ctLP 
= 0.081 and l^d = 0.194. (b) G(1 /bwt 5 Id^p) through the dot at zero l^d- Peaks in 
G correspond to Coulomb blockade on the dot being lifted; the splitting implies finite 
Udg- For hxed Vbwt the difference in peak positions gives the dot-grain charging energy 
Udg = (e)(aLp)(Al/Lp) = 0.081 * 0.26 meV = 21 /reV. Dot-grain tunneling is negligible in 
this limit, (c) Conductance through the grain appears where expected given the interpre¬ 
tation of (b). The conductance is measured with gates BL and BR depleted, measuring 
through the two point contacts formed by gate pairs LBB / BWB and LBT / BWT. 
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adjacent to the gates. The red shaded area in Fig. S3a is 1.2 /rm^. Approximating the 
capacitance of the grain as that of a flat disk with radius r, C = 8er where e = 13eo for 
GaAs and the effective r = 0.62 /rm. This gives an expected Eq = 280 /reV, which is within 
a factor of two of the measurement. In a previous cooldown of the same device we measured 
Ec = 150 /ieV, which we use as the model parameter in NRG calculations. 

In designing the device we aimed for as large an Ec as possible while still being able 
to imagine a near continuum of states in the grain. The level spacing may be estimated 
by considering a particle in a 2D box. The level spacing A = h^7r^/2mA, where A is the 
area of the box and m = 0.067me is the effective mass in GaAs. Using the design area 
A = 1.2 yum^ we hnd A = 4.6 /reV = 2.6/cTe, where = 20 mK. If we instead take 
A = 0.93 yum^ inferred from measurement of Ec and the approximation for the capacitance, 
we hnd A = 6.0 yueV= 3.4/cTe. In either case, A is no more than factors of a few times 
Tg, keeping in mind that the width of the Fermi-Dirac distribution is approximately SAkT^. 
This implies that the grain is indeed acting as a metallic grain at all measured temperatures. 
In Fig. S3, it appears that the typical level spacing (spacing in Ud between diagonal lines) is 
larger than anticipated. The peak conductance can differ signihcantly for each level, which 
rehects a distribution in source-drain coupling asymmetry from level to level. Some levels 
may not be visible if their source-drain coupling asymmetry is strong. 

Measurements of G(Vbwt, Ulp) yield Udg = 21 yueV (Fig. S4). This analysis considers 
the dot-grain system as a capacitively-coupled double quantum dot. The Udg we extract 
should be thought of as an upper bound—the gate voltages are set to a very different regime 
where Fg is negligible, unlike in the paper. When tuning between this regime and the regime 
where Fq ~ F, it appears as if the splitting of the lines in Fig. S4 goes to zero long before 
Fg becomes a signihcant fraction of F, perhaps implying that Udg —t 0. 

S5 Scaling along the quantum critical lines 

In Fig. 2e we demonstrate that measured G{Vsd,T) falls onto a scaling curve derived from 
the conformal held theory (GFT) results of Affleck and Ludwig. In Fig. S5 we demonstrate 
the scaling at other points along the 2GK lines. In most examples, the 2GK scaling behavior 
is faithfully reproduced by the data except perhaps at T = 20 mK. 

A priori, the deviations could be attributed to finite T*, perhaps because data were not 
taken finely enough. However, it seems that the deviations typically appear most pronounced 
on the positive-Ud side. Another possibility could be that true zero Ud drifted slightly over 




Figure S5: Top-left: Measured I^wt) of Fig. 2c. Panels at right: Measured G(14d, T) 

at six points on 2CK lines in the (f^p, Vbwt) plane. Black lines are hts to thermally broad¬ 
ened spectral functions from Affleck and Ludwig (Ref.^°) with small phase shifts from 
potential scattering. 


the course of the measurement. Typically some small applied 14d is required to compensate 
for an offset voltage at the current ampliher input, but at base temperature the quality of 
the scaling collapse is sensitive to errors of just 1 /iV in identihcation of true zero 14d- 


S6 Sensitivity of T* and Sp to fitting range 

In Fig. 3 we use the crossover CFT to £t experimental data and thereby extract the Fermi 
liquid scale T* and the scattering phase shift 6p. The htting procedure uses a limited range 
of 14d (±6 /iV) and it is argued that this is a conservative approach. 

In Fig. S6 we show that the htting is insensitive to small changes in the htting range. At 
each value of Vbwt, we try and extract T* and 6p for nine diherent ranges of bias voltage, 
which we obtain by starting with (-6, -|-6 /iV) and adding or subtracting a point on either 
end, e.g.: (-7.5, -|-6), (-6, -|-6), (-6, -1-7.5), (-4.5, -|-6), etc. It is important not to add so many 
points that data outside the validity of the theory are included. However, subtracting too 
many points may degrade the ht quality. 

After hnding T* and the error in T* reported by the hts, we consider all nine htting 
ranges to give independent estimates of T*, and hnd the weighted mean T*, weighted by the 
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Figure S6: Sensitivity of T* and bp to fitting range. Top panel: source-drain bias spec¬ 
troscopy at T = 20 mK exactly as in Fig. 3. Middle and bottom panels: The black points 
and red curves in the T* and bp panels are exactly as in Fig. 3. The blue points correspond 
to the weighted mean of extracted T* and bp for an ensemble of htting ranges. The error 
bars on the blue points correspond to the standard deviation of the weighted mean. 
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errors from each fit. The error bars show the standard deviation of the weighted mean, and 
indicate the spread of T* values returned by the hts. We do the same for 5p. Varying the 
htting range by small amounts does not seem to contribute signihcantly to the uncertainty 
in T* and 5p. 
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